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FOREWORD 


This  report  was  prepared  by  Dr.  Helmut  Moritz,  Research  Associate, 
of  the  Institute  of  Geodesy,  Photogrammetry  and  Cartography  of  The  Ohio 
State  University,  under  Air  Force  Contract  No.  AF  19(628)^1628,  OSURF 
Project  No.  1474,  under  the  supervision  of  Dr.  Weikko  A.  Heiskanen, 
Director  of  the  Institute.  The  contract  covering  this  research  is  adminis¬ 
tered  by  the  Air  Force  Cambridge  Research  Laboratories,  Office  of 
Aerospace  Research,  Laurence  G.  Hanscom  Field,  Bedford,  Massachusetts, 
with  Mr.  Owen  W.  Williams  and  Mr.  Bela  Szabo,  Project  Engineers. 


ABSTRACT 


General  formulas  for  error  propagation  in  linear 
integral  transformations  are  given.  They  are  applied 
to  the  accuracy  of  the  computation  of  the  gravity  vector 
at  high  altitudes  from  surface  gravity  anomalies,  and 
to  the  accuracy  of  satellite  orbits  computed  therefrom. 


ON  THE  THEORY  OF  ERRORS  OF  PHYSICAL  GEODESY 


(GRAVITATIONAL  FIELD  OF  THE  EARTH  AND  SATELLITE  ORBITS) 


1.  Introduction 

The  computations  in  gravimetric  geodesy  are  mainly  evaluations  of 
some  integrals;  for  computing  satellite  orbits  one  has  to  solve  differ¬ 
ential  equations,  which  again  can  be  reduced  to  integrations.  To  these 
cases,  however,  the  usual  theory  of  error  propagation  cannot  be  applied. 
We  shall  make  this  clear  by  the  example  of  Stokes'  formula,  it 

reads 


f  It 

J  JS(yATA)  4g(9,X)coi9  dy  dX., 

**$—1  X-o 


where 


cp,X  =  geographical  coordinates, 

Q  =  geoid  undulation, 

Ag  =  gravity  anomaly, 

R  =  mean  radius  of  the  earth, 

7  *  mean  theoretical  gravity; 

S(<p,\,9,\)  is  the  function  of  Stokes. 

This  formula  can  be  considered  as  a  linear  integral  transformation 
riich  correlates  the  function  C(cp,\)  to  the  function  Ag(cp,\)  —  both 
gravity  anomaly  and  geoid  undulation  are  functions  of  position  on  the 
sphere . 
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For  deriving  the  accuracy  of  £  from  the  given  accuracy  of  Ag, 
we  have  to  investigate  the  error  propagation  with  respect  to  linear 
integral  transformations.  This  case  is  not  dealt  with  by  the  usual 
theory  of  errors.  To  be  sure,  it  would  be  possible  to  approximate  the 
integral  by  a  sum,  but  this  would  be  unsatisfactory  theoretically  and 
very  complicated  practically. 

In  [4]  or  [5],  however,  we  have  given  a  general  theory  of  errors 
which  can  be  applied  immediately  to  our  case.  It  is  based  on  the  close 
analogy  between  functions  and  vectors  on  one  hand,  and  integral  trans¬ 
formations  of  functions  and  linear  transformations  of  vectors  on  the 
other  hand. 

A  few  results  of  this  theory  which  are  essential  for  the  appli¬ 
cations  to  follow,  will  now  be  given  and  explained,  but  not  derived. 

2.  General  formulas 

In  order  to  facilitate  understanding,  we  first  consider  functions 
of  one  variable,  instead  of  functions  of  two  variables  <p,\.  A  linear 
integral  transformation  of  such  a  function  f(x)  is  of  form 

P 

(1)  f*(x)  =  J  a(x,u)f(u)du 

u<x 

(there  correspond  x,u,f,f*,a(x,u)  to  the  quantities  (cp,\),  (<p,\),  Ag, 

C,  S(cp,\,cp,\),  respectively,  in  Stokes'  formula). 

This  linear  integral  transformation  (1)  is  closely  analogous  to 
the  ordinary  linear  transformation, 


2 


The  error 


(2)  f*  ■  £  .lrfr, 
r=l 


of  the  quantities  which  can  be  considered  as  a  vector 
propagation  for  (2)  is  given  by  the  well-known  formula 


<3)  ■’*“  "A  A v  ^  o«  - 


where  and  (a*^)  are  the  error  covariance  matrices  which  charac¬ 
terize  the  accuracy  of  the  f^'s  and  f* ’ s,  respectively *  E.g,, 


ik 


Mi 

013  • 

»  •  °in\ 

°22  • 

•  •  a2n| 

Wi 

0X12  «  ' 

■  •  Jnn  / 

The  diagonal  terms  are  nothing  other  than  the  squares  of  the  standard 
errors  of  f^j 


(4) 


while  0  ,  i^J,  expresses  the  mutual  correlation  of  f  and  f  ---  if 

i  J  i  j 

all  f^  are  independent,  all  cr^ ^ )  are  zero. 


Cf *  Grossmann  [l],  pp.  333  and  334;  his  formulas  (21)  and  (23)  are 
equivalent  to  our  eqs„  (2)  and  (3),  since  the  matrix  of  cofactors,  Q, 
is  proportional  to  the  error  covariance  matrix 0  ^o^ik*  wtiere 

is  the  standard  error  of  unit  weight. 
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1 


i 

I 


(2)  and  (1)  are  completely  analogous,  there  correspond 


function  f(x)  to  vector  f  j 
(5)  variables  x,u  to  indices  i,r; 

P  n 

integral  J  to  sum  £  . 

u<x  r=l 


Hence  it  is  evident  that,  corresponding  to  the  error  covariance 
matrix  which  is  a  quantity  with  two  indices,  we  shall  have  a  function 

»  i 

a(u,  u  )  of  two  variables  u,u  ,  which  will  be  called  error  covariance 
function,  or  error  function.  Corresponding  to  (4)  we  have 


(6)  a(u,u)  =  m3(u), 

»  I 

the  square  of  the  standard  error  m(u)  of  f(u),  while  a(u,u  )  for 

1 

characterizes  the  mutual  correlation  of  the  values  f(u)  and  f(u  )t 

Since  a(u,u  ),  in  most  cases,  is  a  continuous  function,  it  is  even 
for  u^ur,  in  general,  different  from  zero,  sc  that  there  is  always 
correlation*  It  is,  however,  often  possible  to  assume  m(u)  to  be  inde¬ 
pendent  of  u  and  constant  *  In  this  case,  cr(u,u')  is  a  function  of  the 
difference  u-u*  only*  This  function  often  has  a  shape  similar  to  Fig*  1 
only  the  neighboring  values  are  correlated,  the  correlation  being  a 
decreasing  function  of  distance  |  u-uf|  « 


£(u,u') 


Fig  *  1 
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Functions  similar  to  Fig.  1  are,  e.g., 


(7.)  o<«,u’>. 

or 

»  ° 

(7b)  o(u,u  )=  i;+ka(u!u<)a  • 

From  (3),  t>y  means  of  the  analogies  (5),  we  find  the  formula  of 
error  propagation  in  linear  integral  transformations  (l): 

i  &  it  i  i 

(8)  o*(x,x  )  =  j  j  a(x,u)  a(x  ,u  )  c(u,u  )  du  du  . 

u=a  u’=a 

We  get  an  important  special  case  if  we  need  the  function  f*(x)  in 
one  fixed  point  x=xq  only.  Here  we  introduce  the  following  notations 

f*(xQ)  s  F,  o*(xQ,x0)  =  m*a(xo)  a  M3,  a(xQ,u)  •  h(u)  . 

F  Is  a  linear  functional  of  f(x);  according  to  (1)  it  is  given  by 

e 

(9)  F  =  J  h(u)  f(u)  du  . 

u=a 


The  corresponding  formula  of  error  propagation  is 


(11) 


Using  (5),  we  have  in  the  analogous  case  of  the  vector  f  - 


n 


and 

n  n 

(12)  M®=  E  Eh  h  a  , 

x  '  .  -  r  s  rs 

r=l  s=l 

i.e.,  the  well-known  elementary  case  of  a  linear  function  of  the 
variables  fj,...f  and  its  error  propagation*  So  our  linear  functional 
simply  corresponds  to  a  linear  function. 

Between  (1)  and  (2)  and  between  (9)  and  (ll)  we  have  further 
analogies  which  are  still  closer.  If  we  develop  a  function  f(x)  in 
a  series  with  respect  to  a  system  of  orthogonal  functions,  the  coef¬ 
ficients  of  this  development  form  an  infinite  vector;  and  if  we  per¬ 
form  the  same  development  for  the  error  function  a(x,x*),  the  coef¬ 
ficients  form  an  "infinite  matrix."  With  these  infinite  quantities 
we  can  operate  in  exactly  the  same  way  as  with  finite  vectors  and 
matrices.  Hence,  eqs.  (2),  (3),  (11),  and  (12)  hold  immediately  for 
this  case;  we  only  have  to  replace,  in  the  summation  sign,  n  by  ». 

Thus,  any  linear  integral  trans format ion  is  equivalent  to  a  linear 
transformation  of  the  coefficients  of  an  orthogonal  expansion. 

Such  orthogonal  expansions  are,  e.g.,  Fourier  series  or  expansions 
in  spherical  harmonics.  They  are  important  because  a  very  simple  trans¬ 
formation  of  the  coefficients  may  correspond  to  a  complicated  integral 
transformation.  A  good  example  is  Stokes*  formula;  this  fact  is  used 
in  deriving  this  formula  by  means  of  spherical  harmonics. 

Since  spherical  harmonics  are  of  fundamental  importance  in  phy¬ 
sical  geodesy,  the  error  theory  of  orthogonal  functions  is  Important, 
too.  In  the  examples  treated  in  this  paper  we  do  not  need  them,  how¬ 
ever.  Therefore  we  cannot  present  them  here.  The  reader  will  find  a 
comprehensive  treatment  of  the  whole  subject  in  [51* 

3.  Accuracy  of  the  Gravity  Vector  at  High  Altitudes 

Let  us  now  apply  these  general  formulas  to  problems  of  physical 
geodesy.  We  assume  gravity  on  the  earth’s  surface  given  with  a  certain 
accuracy.  From  these  data  we  compute  the  gravity  vector  outside  the 
earth.  By  means  of  this  external  gravity  vector  we  can  calculate 
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satellite  orbits.  How  accurate  is  the  computation  of  the  external 


gravity  vector  and  of  satellite  orbits?  We  shall  simplify  our  assump¬ 
tions,  in  order  to  make  clear  the  structure  of  the  problem. 

Besides  the  actual  gravity  vector  ,  we  introduce  the  theoret¬ 
ical  gravity  vector  ^  ,  referred  to  the  same  point,  and  the  vector 

of  gravity  disturbance,  5  ,  so  that 

-* 

and 

(13)  S  =  -grad.  T  . 

Hence,  5  is  the  gradient  vector  of  the  disturbing  potential  T  which  is 
defined  as  the  difference  of  actual  gravity  potential  W  and  theoretical 
gravity  potential  U. 

T  may  be  expressed,  according  to  Helmert,  by  a  "coating”  with 
density 

(1*0  d  =  Ag  +  ||  Q, 
namely, 

(»>  T  -  5,  XT  5  “• 

The  integration  is  performed  overthe  globe;  D  is  the  distance  of  the 
surface  element  dS  from  the  point  P  in  which  T  is  computed. 

In  our  case  it  is  sufficient  to  replace  the  globe  by  its  tangential 
plane,  so  that  by  (15), 
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I 


i 

*: 

(15’)  T(x,y)  =  ~  6JW  JL  d*dn  • 

Here  is 

(16)  D3  =  (x-|)3  +  (y-r,)3  +  H3  ; 

x,y  are  plane  rectangular  coordinates  of  point  P  ;  H  is  the  height 
above  ground;  and  g,  r\  are  the  coordinates  of  dS. 

Denoting  the  components  of  the  vector  of  gravity  disturbance  by 

v  v  v we  find  by 

,  _  dT  ,  dT  ,  _  dT 

x  "  "dx»  V  'dy*  V  "dH  • 

i  I 

Differentiating  (15  )  with  respect  to  x,  y,  H  yields 

-  oo 

(17) 

-oo 

These  formulas  are  integral  transformations  which  correlate  the 

functions  6.6.6  to  the  function  u* 
x’  y7  z 

For  applying  the  general  formulas  of  the  foregoing  section  to  this 

8 


f 


problem,  we  have  to  consider  that  here  we  have  functions  of  two  variables, 

the  error  functions  of  which  are,  therefore,  functions  of  four  variables. 

The  coordinate  pairs  (x,y)  and  correspond  to  the  single  variables 

x  and  u,  in  the  previous  section;  instead  of  a  single  integral  we  now 

have  a  double  integral.  Denoting  the  error  functions  of  6  (x,y),  6  (x,y), 

x  y 

»»  i  i  »i 

6z(x,y)  and  of  n(x,y)  by  ax(x,y,x,y  ),  ay(x,y,x  ,y  ),  az(x,y,x  ,y  ), 

I  I 

and  a(x,y,x  ,y  ),  respectively,  we  find  by  (8)  immediately 


(18)  ~  fjff  dfdr)', 

-  oo  V  J 


drd,,/- 


where 


(7a) 

(19) 


D  is  given  by  (l6)  and  D  by 

d's  =  (x'-g‘)a  +  (y'-T)'  )a  +  H3  . 

For  simplicity,  we  assume  the  error  function 
(two-dimensional) : 


-66e' 


[C^-^/)1+(,?-rl/)2J 


a  to  have  the  form 

2  yt 
-cd 

e 
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i.e.,  the  correlation  depends  only  on  the  distance 


d  =  +  (n-V)3 

t  i 

of  points  ( |, rj)  and  (§  ,r\  )  <,  (19)  being  appreciably  different  from  zero 

t  1 

for  (|,i|)  =  (|  , t]  )  only,  we  may  set  in  (l8),  approximately, 

I  t 

I  =  I,  n  • 

So  we  get,  say,  for  a 


-c2d2 


•&!!!!  W 

-00  4)  ' 


where 


D"a=  (x  -|)3  +  (y  -r])3  +  H3 


Now  we  first  integrate  with  respect  to  |  ,  t)  ,  where  x  -{,  D  and 

D"  are  constant  for  the  integration.  Since 

-Ml 


OO 


1  / 


£'--0O  /)'*-*> 

©o 

*  <5*0  *  J  £ 

-OO 

-  ir-Sfc- 
cz 

we  obtain 

(20)  erx(x,y 

oO 


d*d 
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Even  now  direct  integration  is  very  difficult,  resulting  in  elliptic 
integrals . 

VJe  may,  however,  avail  ourselves  of  a  development  in  a  series.  By 
S-x  ~  cos  cp,  rj-y  -  r  sin  cp 

and 

;21)  x  -x  -  s  cos  a>  y  “y  =  s  sin  a 


ve  introduce  polar  coordinates,  the  origin  being  point  P(x, y)  (Fig.  2). 


->  X 


Fig.  2 

Then, 

D*'  —  H2-r2,  B"a=  H8+i,,3=  Fa~r2-  2rs  cos(cp-a)  +  s2  . 
Developing,  for  small  s,  m  a  power  series  with  respect  to  s,  we  get 

i  |  ,  3 r  ccn  (<p -oO  ,  ;  3  15  r2  coF(<f-oC)l  cC 

F5  '  F  +  — P — 5  l'  IF  T  X - 5^— F  +  • ' 
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Inserting  this  in  (20)  and  considering 


x- £  s  -r  cos  cp, 


x  -|  =  -r  cos  cp  +  s  cos  a  , 


b  =  rr^)3/2 ' 


d|dr)  s  r  dr  dcp  , 


we  may  perform  the  integration  in  polar  coordinates  relatively  easily. 
In  a  similar  way  we  treat  and  a finally  obtaining 


Oo 

r,  3  s2 

16c*HT 

-1  16  TF 

1 

1  ±  s2 

16caHM 

.  16  TF 

6ro  [■ 1 

3  s2 

8  TF  ■ 

•j. 


The  polar  coordinates  s,  a  are  related  to  x,  y,  x  ,  y  by  (21). 

We  see  that  the  expressions  (22)  depend  on  the  mutual  position  of  points 

i 

P  and  P  only:  a  depends  only  on  the  distance  s,  a  and  a  on 

^  y 

the  direction  a,  too. 

According  to  (6)  we  get  the  mean  square  errors  of  6  ,  6  ,  6  by 

x  y  z 

i  t 

setting  in  (22),  x  =  x,  y  =  y,  i.e.,  s=0.  Hence, 


(23)  m3  «  m3  =  - - — ,  m3  =  - — 

x  y  Z  Sc3!!3 


mx=m^  and  have  the  same  order  of  magnitude;  they  are  inversely 
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proportional  to  the  elevation  H  above  ground.  These  formulas  become 
invalid  for  small  H  because,  in  this  case,  (23)  tends  to  infinity. 
For  small  H,  indeed,  the  simplifications  leading  to  (20)  would  not 
have  been  admissible. 

More  detailed  investigations  show  that  the  error  functions  (22) 
can,  in  sufficient  approximation,  be  represented  by  functions 

I+k'P'  =  m3(l-kas3  +  ...)* 

Thus  we  find  from  (22), 


1 

mx 


) 


(24) 


^(S,OC) 


mv 


1  +■  ^  ( 1  +■  2  M/n.  oc) 


er2(s)  - 


m. 


1  + 

1  8  TT 


These  are  the  desired  error  functions  of  6  ,  6  ,  6  ,  or  of  the  eompo * 

x  y  z 

nents  gy,  gz  of  the  gravity  vector  -oy  . 
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4.  Influence  on  Satellite  Orbits 

In  the  last  section  we  have  evaluated  the  accuracy  of  the  gravity 
vector  (24);  now  we  are  going  to  estimate  the  influence  of  this  accuracy 
on  computed  satellite  orbits „ 

Essential  for  the  following  are  the  formulas  for  the  change  of 
the  orbital  parameters  a,  e,  i,  Q,  caused  by  a  disturbing  force 
(X,  Y,  Z),  which  are  found  in  any  textbook  on  celestial  mechanics 


(cf.  [8]): 


da 

dt 

de 

dt 

di 

dt 

dfl 

dt 


2a  h[X  sec  cp(l  +  e  cos  v)  +  Z  tg  cp  sin  v], 


h  cos  cp[x( cos  v  +  cos  E)  +  Z  sin  v], 


r  h 


a  cos  cp 


Y  cos  u, 


r  h  y  sin  u 
a  cos  cp  sin  i 


Here  we  have  used  the  following  notations: 

a  s  semi -major  axis  of  the  satellite  orbit, 
e  =  numerical  excentricity  of  the  orbit, 
i  =  inclination  of  the  orbit, 

Q  =b  right  ascension  of  the  ascending  node, 
t  =  time 

r  =  radius  vector 
v  =s  true  anomaly 
E  =  excentric  anomaly 

u  =  longitude  in  the  orbit,  reckoned  from  the  node  (cf0  [3])» 
cp  is  an  auxiliary  angle  connected  with  the  excentricity  e  by  the 
relation 

e  =  sin  cp  . 

h  is  a  constant  defined  by 
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{  f  =  Newton's  constant  of*  gravitation,  M  =  mass  of  the  earth).  X,  Y,  Z 
are  the  components  of  the  disturbing  force,  namely 

Z  directed  along  the  radius  vector, 

Y  normal  to  the  orbital  plane, 

X  normal  to  Y  and  Z  , 

We  restrict  our  problem  to  the  simplest  and  practically  important 
case  of  an  approximately  circular  orbit.  Here, 


e  =  sin  cp  ±  0,  r  ^  a,  E  *  v, 


and  we  have 


=  2a  h  X, 
at 


(25) 


de 

dt 


=  2  kX  cos  v  +  h  Z  sin  v, 


di 

at  =  H  Y  cos  u> 


4Q  „  sin  u 

dt  =  kYhh-i 


X  is  now  directed  along  the  tangent  to  the  orbit;  X  and  Y  are 
horizontal,  Z  is  vertical.  Hence,  we  can  identify  the  directions  of 
the  coordinate  axes  x,  y,  z  of  the  foregoing  section,  with  the  directions 
of  X,  Y,  Z.  Then,  the  components  X,  Y,  Z  of  the  disturbing  force  are 
nothing  other  than  the  gravity  disturbances  6^,  6^,  6^. 

We  consider  first  da/dt.  Denoting  by  Aa  the  finite  variation  of 
a  for  a  whole  revolution  (period  of  revolution  =  T),  we  obtain 


2jr 


t 


in  ; 


t=0 


u=0 


u=0 


da 

dt 


du, 
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hence 

t 

Aa  b  2a  *  --  J  X(u)  du  . 

^  0 

This  is  a  linear  functional  of  form  (9)(f(u)  =  X(u),  h(u)  = 
2a  h  T/&r  =  const.)*  Eq.  (10)  yields  immediately 

_s  2tt  2tt  ,  , 

(26a)  M3  =  4a3n3  -  f  fa  (u,u  )  du  du  . 

a  4tt3  u=0  u'=0  x 


Analogous ly, 


rf  Zrr 


m3  =  ^  a  P  fa  ( v, v  )  cos  v  cos  v  dv  dv  + 

*  ^viovdo  x 


(26b) 


+  K 


rpS  ^  » 


J  J  CT  (v,V  )  Sin  V  Sin  V 


f  » 


Ett3’  o  z 

H7r  v=0  v’=0 


dv 


(since  6  and  6  can  be  shown  to  be  practically  unc or related) ,  and 
x  y 

a  T3  *  •  i 

(26c)  Mr  *  h  t— ■ t  J  J  a  (u,u  )  cos  u  cos  u  du  du  , 

1  T  u=0  u'=0  y 


2/r  2jr 


m*  *£«  I  II 

(26d)  M?  =  -T-av  rar  J  J  a  (u,u  )  sin  u  sin  u  du  du 
w  sin  1  47r  u=0  uf=0  y 


The  error  functions  a  ,  az  are  given  by  (24).  The  x-axis 
being  tangent  to  the  orbit,  we  have  along  the  orbit 
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a  =  0,  s  =  a(u-u  )  =  a(v-v  ) 
(u  and  v  axe  angles).  Hence, 


CTx(u'u  }  =  IsA3-  [1  +  iIf  (u'u  )3]  1  ’ 

CTy(u,u  >  =  I&V  [1  +  i§F  (u"u  )3]  1  ’ 


oz(u,u  )  = 


BS«1F  [1  +  85*  (u-u  )3]  1  • 


Inserting  this  in  eqs.  (26)  and  applying  a  similar  trick  as  in 
integrating  (18),  we  find  in  sufficient  approximation  finally 


(27) 


<  =  i 


w3  _  ^+/6  a 

“fe  =  ~W  Q  ’ 


/3 


3  _  0“ 

i  _  CB  y  ’ 


vt 


1  /3  03 

‘Q  slnai  58  y  ’ 


where 

Qa  =  —  SO  l 
W  tM  ?H 


is  a  nondimensional  quantity. 


By  (27)  we  have  immediately 


,  M  =  0.408  aQ,  =  O.367  Q, 

(27  )  a 

M  =  0.190  Q,  =  0.190  Q/sin  i  . 


These  are  the  mean  square  errors  of  the  variations  Aa,  Ae,  Ai,  AQ 


for  one  revolution,  which  are  caused  by  the  inaccuracies  of  the  given 


gravity  values,  measured  on  the  earth* s  surface.  They  are  valid,  provided 


17 


that  the  error  function  of  the  coating  m*  can  be  represented  with 
sufficient  accuracy  by  a  function  (19)  with  parameters  gq,  c. 

In  order  to  get  the  mean  square  errors  of  the  variations  of  orbital 
elements  for  several  revolutions,  we  obviously  have  only  to  multiply  the 

1 

values  (27  )  by  the  square  root  of  the  number  of  revolutions, 

5.  A  Numerical  Example 

We  have  computed  some  numerical  values  for  an  elevation 
H  =  500km 

above  ground.  For  and  c  we  choose  some  characteristic  pairs  of 
values.  The  results  are  shown  in  Table  1. 


Table  1 


a 

0 

mgal 

c-1 

km 

m  =  m 

Wi  y 

m 

mg&l 

Ma 

m 

Me 

Mi  =  ^  Sin  1 
seconds  of  arc 

2 

15 

0,01 

0.02 

0.3 

0.04x10”^ 

0.005 

10 

25 

0.1 

0.2 

3 

O.UxlO-6 

o.ou 

15 

45 

0.3 

0.5 

7 

l.OxlO-6 

0.11 

19 

90 

0.9 

1.2 

18 

£ 

2.5x10 

0.27 

20 

220 

2.2 

3-1 

k6 

6.6xl0‘6  ■ 

O.69 

The  pairs  /o  ,  c"  have  the  following  meaning.  We  make  the  simpli¬ 
fying,  but  not  qu?te  correct,  assumption  that  the  inaccuracies  in  M-(l^) 
are  caused  only  by  Ag,  but  not  by  £.  Then,  interpolation  being  the 
main  source  of  error  in  Ag,  a  is  the  error  function  for  interpolation 
of  gravity,  which  is,  in  fact,  of  form  (19)'  /crQ  is  the  mean  square 
error  of  integration,  and  c  has  the  dimension  of  a  reciprocal  length, 
so  that  /aQ  and  c*1  are  more  intuitive  than  aQ  and  c.  The  pairs  of  values 
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given  may  correspond  to  the  interpolation  error  for  average  station 
distances  of  about  10,  50,  100,  200  and  500  km  [7]** 

The  values  found  in  this  way  are  very  small.  Even  if  we  correctly 
have  regard  for  the  errors  of  the  geoid  undulations,  its  order  of 
magnitude  will  not  be  changed,  as  can  be  seen  on  closer  examination. 

A  necessary  condition  is,  however,  uniform  coverage  of  the  whole  earth 
by  gravity  stations. 

We  have  deliberately  simplified  conditions  as  much  as  possible 
because  of  the  difficulty  of  the  problem.  More  realistic  estimations 
will  have  to  take  into  account,  e.g.,  the  errors  of  Q  and,  above  all, 
the  big  gaps  in  the  world  gravity  net  which  considerably  impair  the 
accuracy. 


*  These  values 
by  Kaula  (J. 


are  consistent  also  with  the 
Geophys.  Res.  64,  2401-2421, 


covariance  function  given 

1959) • 
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